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ABSTRACT 

We use multi-epoch spectroscopy of ^4000 white dwarfs in the Sloan Digital Sky Survey to constrain 
the properties of the Galactic population of binary white dwarf systems and calculate their merger 
rate. With a Monte Carlo code, we model the distribution of ARVmax, the maximum radial velocity 
shift between exposures of the same star, as a function of the binary fraction within 0.05 AU, /bin 
and the power-law index in the separation distribution at the end of the common envelope phase, 
a. Although there is some degeneracy between /bin and a, the the fifteen high ARVmax systems 
that we find constrain the combination of these parameters, which determines a white dwarf merger 
rate per unit stellar mass of 1.4^^ q x 10""'^'^ yr~^MQ^ {la limits). This is remarkably similar to the 
measured rate of Type la supernovae per unit stellar mass in Milky Way-like Sbc galaxies. The rate 
of super-Chandrasekhar mergers is only LO^J g x 10~^^ yr~-'^MQ^. We conclude that there are not 
enough close binary white dwarf systems to reproduce the observed Type la SN rate in the 'classic' 
double degenerate super-Chandrasekhar scenario. On the other hand, if sub-Chandrasekhar mergers 
can lead to Type la SNe, as has been recently suggested by some studies, they could make a major 
contribution to the overall Type la SN rate. Although unlikely, we cannot rule out contamination 
of our sample by M-dwarf binaries or non-Gaussian errors. These issues will be clarified in the near 
future by completing the follow-up of all 15 high ARVmax systems. 

Subject headings: binariesxlose, spectroscopic — white dwarfs — supernovae: general 



1. INTRODUCTION 

The nature of the progenitor systems of Type la su- 
pernovae (SN la) remains one of the key open issues in 
stellar evolution. There is a general agreement that the 
exploding star is a CO white dwarf (WD) that is some- 
how ignited following accretion of material from a bi- 
nary companion, but the identity of this companion is 
still a matter of debate. Most theoretical scenarios for 
SN la progenitors can be divided into two broad classes : 
single degenerate (SD) systems (jWhelan fc IbenI flQTSf ). 
where the companion is a non-degenerate star, and ma- 
terial is transferred over < 1 Myr via Roche lobe overflow 
or wind accretion and double de gener ate (DD) systems 
(|Iben fc Tutukovl 119841 : IWebbinklflQSl . where the com- 
panion is another WD, and mass transfer happens over 
much shorter time scales in a merger event. 

Recent developments have provided some evidence 
in favor of the DD scenario. Sensitive searches have 
failed to find signs of a mass-losing nondegener- 
ate co mpanion in ra di o ob servations (jHoresh et al.l 
120111: iChomiuk et all I2012D. early l ight curves 
(IHavden et al.l 120101: IBianco et alj 120111: Bloom et al.l 



2007 



20111) and supernova remn ants (iBadeiies et al . 

Schaefer fc Pag notta '2012' but d. I Williams et"!!] 
and see [Stcriibcrg e t all 120111 for some statistical 



evidence of absorption in SN la spectra that could be 
of SD origin). Several independent measurements of 
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the delay time distribution (the specific rate of SN la 
as a function of time after a hypothetical brief burst 
of star formation) have converged onto a. ^ shape 
extending from a few hu ndred Myr to several Gyr (see 
iMaoz fc Mannuccil [20TTl for a recent review), a form 
expected from a population of DD systems that merge 
due to gravitational wave emission. 

In this Letter, we present the first measurement of the 
local WD merger rate and use it to test the viability of 
the DD scenario for SN la progenitors. To date, orbital 
parameters have been measured for over 40 individual 
DD systems with p eriods ranging betwee n 12 minutes 
(iBrown et al.|[2011bD and several days (e.g. iMarsh et al.l 
119951 iNelemans et al.l [2005t iKilic et al.ll2010D . but the 
fundamental properties of the Galactic DD population 
are still poorly known. Some small WD samples have 
been used to estima te parameters like the binary frac- 
tion (jMaxted fc Mar sh 1999, find it to be between 0.017 
and 0.19, based on 46 WDs) or the merger rate limited 
to systems that contain extremely low mass (< 0.25 Mq) 
WDs ( ^ 4 X 10~^ yr~^ based on 12 svstems. lBrown et al.l 
l2011al) . Here, we adopt a statistical approach to char- 
acterize the DD population, using the unique capabil- 
ities of the Sloan Digital Sky Survey (SDSS) for time 
resolved spectroscopy of thousands of WDs. We use this 
large data set to measure the distribution of ARVmax, the 
maximum radial velocity (RV) shift between exposures of 
the same WD. In a companion paper (Maoz, Badencs fc 
Bickerton 2012, henceforth Paper I), we describe how this 
distribution constrains the binary fraction and the sep- 
aration distribution. Because the evolution of detached 
DD systems is driven only by gravitational wave emis- 
sion these parameters, together with the WD masses, 
uniquely determine the local merger rate, which can be 
compared to measurements of the specific SN la rate in 
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nearby Milky Way-like galaxies. 

2. OBSERVATIONS 

The SDSS (|York et al.ll2000[ ) co ntains the largest avail- 
able collection of W D spectra (jKleinman et al.l 120041 : 
lEisenstein et al.ll2006[) . The latest version of the SDSS 
WD catalog, corresponding t o DR7, has over 17000 en- 
tries emman et al.ll200m For simplicity, we have re- 
stricted our analysis to non-magnetic DA WDs that show 
no obvious signs of main sequence companions, and have 
no confirmed or possible absorption lines from any ele- 
ments other than H - i.e., objects classified as 'DA' or 
'DA:' in the catalog, We further removed from this sam- 
ple seven WDs with main sequen ce companions listed by 
iRebassa-Mansergas et al.l (|2010( ). leaving 12763 objects. 

Like all SDSS spectra, the spectra of these WDs were 
divided into three o r more sub-exposures to facilitate cos- 
mic ray rejection (iStoughton et al.ll2002[ ). Since DR7 
(|Abazaiian et al.l I2009D . these sub-exposures are avail- 
able from the SDSS server. The ongoing SWARMS 
survey is using these sub-exposures to identify short- 
period binar y WD candidates in SDSS, which are the n 
followed up (|Badenes et all[2009t iMuUally et al.l 120091). 
The DS/DT collaboration (jBickerton et al.ll20nT has im- 
plemented a pipeline to handle cosmic ray rejection and 
derive stable wavelength solutions for all sub-exposures 
of the same object. This Letter presents the first results 
obtain ed using the DS/DT pipeline - see feickerton et al.l 
(|2011| ) for technical details. 

To measure RVs, we normalize each spectrum, divid- 
ing it by a highly smoothed version of itself, and we fit 
four Balmer lines - Ha through - with Voigt pro- 
files in absorption. We find the best-fi t model for each 
line using MPFIT (|Markwardtl l2009( ). an IDL imple- 
mentation of the Levenberg-Marquardt algorithm. For 
each sub-exposure n, these fits yield four RV measure- 
ments RV„,i, and four la errors 6¥(Vn,i [i = 0..3). 
To define a single RV for each sub-exposure, we im- 
pose a threshold on the error of individual RV measure- 
ments, (JRVthr, and discard values with larger errors or 
that deviate by > la from the weighted mean, RV„ = 

{Y.i^n,i X ^n,i) /Y.i^n,i, with Wn,i = l/<5RV„,i. In 

sub-exposures with at least two surviving RVs, we re- 
calculate the weighted mean RV„, and define its error as 

(SRV„ = \Jl/^i w'^ j. Finally, we define ARVmax as the 

difference between the highest and lowest RV„ in each ob- 
ject. This procedure eliminates noisy sub-exposures with 
internally inconsistent RVs, effectively culling most WDs 
fainter than g ~ 19. We find that (5RVthr=80 kms^^ 
is a good compromise between sample size and quality, 
with 4063 WDs and 15,236 RV measurements. The dis- 
tribution of RVs and (5RVs in this sample is shown in 
Figure m The mean of the RV distribution, 31.0 ± 0.4 
kms~^, is non-zero and positive due to gravitational red- 
shift at the WD surfaces. This is the expected value for a 
local WD population with an average mass of ~ 0.6 Mq, 
and confirms the 3 2.6 ± 1.2 kms~^ redshift obtained by 
iFalcon et all (|2010i ) using a different sample of 449 DA 
WDs observed at high resolution (^16 kms~^ per pixel. 

The pre-publication catalog, July 2010 version, was kindly pro- 
vided by S. Kleinman (private communication), and has 17,371 
entries. 



see iNapiwotzki et al.ll200ll - for comparison, SDSS spec- 
tra are 70 kms~^ per pixel). The agreement shows 
that our procedure yields well-calibrated RVs. 

A key ingredient in the Monte Carlo models described 
in Section[3]and Paper I is the distribution of errors, (5RV. 
We therefore estimate the errors independently, and com- 
pare them to the errors reported by MPFIT. For each 
Balmer line of each WD, we find the difference between 
random pairs of RVs that have similar reported values of 
(5RV. These differences have an approximately Gaussian 
distribution, with a a that gives an empirical estimate 
of the true error in the RV measurements, times V2. 
Using this procedure, we find that MPFIT tends to un- 
derestimate the RV errors. We obtain the corrected (5RV 
distribution shown in Figure [1] by running a grid of such 
pair subtraction tests. 

3. THE WD MERGER RATE 

The distribution of ARVmax values in the (5RVthr=80 
kms~^ sample is shown in Figure [5] The vast major- 
ity of WDs have low ARVmax - these are either iso- 
lated WDs or DD systems where the data show no ev- 
idence for significant RV shifts, because the period is 
too long, or the inclination is too small, or the sub- 
exposures were taken at similar orbital phases. However, 
there is a significant tail of WDs with high ARVmax, ex- 
tending beyond 500 kms~^. From the distribution of 
SDSS sampling times and RV errors, we do not expect 
any non-binary WDs to have ARVmax higher than ^300 
kms^^ (dashed line in Figure [2] - see Paper I for de- 
tails). The eight objects with higher ARVmax are there- 
fore real binaries, detected at high significance. Some of 
these are already pu blished discoveries, like SDSS 1257 
dBadcnes et al."2009', ARVmax= 538 kms'^), or SDSS 
0923 (Brown et al. 2010, ARVmax= 410 kms'^). Others 
will be the subject of forthcoming papers. 

At intermediate values of ARVmax, between ^200 and 
300 kms^^ follow-up observations are necessary to tell 
which of the several dozen WDs are real binaries, but 
here we are only concerned with the statistics of the 
ARVmax distribution, not the binary character of indi- 
vidual systems. As a precaution, we have visually vetted 
the spectra of the 107 WDs with ARVmax> 175 kms~^ 
to check for any misclassified objects, WDs with weak 
line emission, spurious line fits, etc., and removed them 
from the distribution. 

The Monte Carlo techniques that we use to model the 
ARVmax distribution are described in detail in Paper I. 
The method involves generating a large number of WD 
systems with a given binary fraction, drawing the bi- 
nary separation, a, and the component masses from given 
distributions, assigning random inclinations and phases, 
choosing the photometric primary, and folding the re- 
sulting RV curves through the observational parameters 
(sampling times and (5RV distribution) of the SDSS data. 
Because we can only detect binaries at high confidence 
with ARVmax^ 250 kms~^, there is an upper limit to 
the separations that we can probe. For an extreme-mass- 
ratio WD pair with Mi = 1.1 Mq and M2 = 0.2 Mg, an 
RV curve with a peak-to-peak amplitude of 250 kms~^ 
corresponds to a separation of ~ 0.05 AU. Therefore, we 
define the measured binary fraction, /bin, as the fraction 
of WDs with companions within Omax = 0.05 AU. As ex- 
plained in Paper I, our definition of /bin does not include 
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Fig. 1. — Left: Distribution of RVs. The nonzero mean (31.0 
Illustration of the pair subtraction test used to correct the errors 
i5RV before (dashed blue line) and after (solid red line) correction. 

an additional population of extremely low-mass binaries 
with primaries below 0.25 Mq, whi ch we assume are al- 
ways in short-period binaries (see iBrown et alj l2011al . 
and references therein). To constrain the distribution 
of binary separations, we model it as a power-law with 
index a at the end of the final common-envelope episode, 
n[a) oc a". In Paper I, we derive analytically the time- 
evolved form of this distribution due to gravitational 
wave emission and merger events, integrated over the 
star-formation history of the Galaxy. This evolved sep- 
aration distribution is used for generating the simulated 
WD pairs. Primary masses arc chosen from t he ob served 
distribution of isolated WDs in lKepler et all ()2007[) . Sec- 
ondary masses are chosen from a power-law mass-ratio 
distribution of index (3. In practice, the results have a 
weak dependence on /3 (see Paper I for details), and the 
analysis we present here assumes (3 = 0. We choose the 
photometric primary by random draw, except for < 0.35 
Mq secondaries, which we assume are always photomet- 
ric primaries (see Paper I). 

Each combination of /bin and a leads to a predicted 
ARVmax distribution in the survey. Examples of such 
model distributions arc plotted alongside the observed 
one in Fig. [2l This plot showcases the discriminating 
power of the ARVmax distribution, in that some combi- 
nations of /bin and a are allowed by the data, and some 
are clearly ruled out. The validity of each model can 
be quantified by multiplying the Poisson probabilities of 
finding the observed number of systems in each ARVmax 
bin, given the expectation from the model. In all com- 
parisons between models and observations, we use only 
the ARVmax bins above 250 kms~^ (15 systems in total). 

The likelihoods of different WD population models arc 
shown in Figure [3l with the 95% and 68% confidence 
level contours marked. Lower binary fractions require 
more negative values of a to reproduce the observations, 
and vice versa. This is because the tail of WDs at high 
ARVmax (i-e. short-period binaries) can result either 
from WD populations where many binaries are formed 
at all separations (high /bin, more positive a), or few 
binaries are formed but with preferentially small sepa- 
rations (low /bin, more negative a). Some population 
models can be ruled out using simple arguments. We can 
require, for instance, that <100% of WDs have a com- 
panion with P<12 days (the circularization limit), which 
woul d be in conflict with m ain sequence binary surveys 
(e.g. iRaghavan et al.|[20Tol ). We can also require that 
WDs do not merge faster than they are formed, because 
we do observe pre-merger systems. These two conditions 
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Fig. 2. — Distribution of ARVmax with Poisson error bars (black 
histogram), compared to model distributions (blue, red, orange 
and green solid curves). The dashed black curve is a model with 

/bin = 0. 

rule out the striped regions in the /bin, ct plane marked 
in Figure [3] 

Every simulated system has aknown merger time, and 
therefore it is straightforward to obtain the merger rate 
for each simulated binary population, whether in total, 
or for specific cases, such as super-Chandrasekhar sys- 
tems. The merger rate per WD in the sample can be 
converted to a merger rate per unit stellar mass in the 
solar neighborhood us ing measured est imates of the local 
WD nu mber densitv fSion et al.ll2006[ ) and stellar mass 
density (|McMillanll2011l ). Curves of constant WD merger 
rate per unit stellar mass are shown in Figure [H As ex- 
plained in Paper I, the curves appear as straight lines in 
a — log(/bin) space. The slope of the solid lines (all merg- 
ers) is l/log(amax/ao) ~ 1-5, where 0.05 AU, 
and oo ~ 0.011 AU is the initial separation of a typ- 
ical WD binary that merges in 10 Gyr (our assumed 
age of the Galaxy; see Paper I). Offsets in a between 
merger rates separated by a decade are also 1.5. Despite 
the broad range of allowed a and /bin values, the SDSS 
data put stronger constraints on the specific WD merger 
rate. The likelihood-weighted merger rate is 1.4 x 10"^"^ 



-14 V1--1 m;:;^ for 



yr Mq for all systems, and 1.0 x 10 yr" i\±q 
super-Chandrasekhar systems. Table [T] gives 95 % con- 
fidence levels on these rates, and allowed ranges of /bin 
for specific values of a. 

4. DISCUSSION AND CONCLUSIONS 
iLi et al.l (|2011f ) have recently measured the SN la rate 
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Fig. 3. — Likelihood contours in the /bim c« plane, indicating 
the 95% (light gray) and 68% (dark gray) confidence levels. The 
overlaid color lines are curves of constant WD merger rate, total 
(solid) and supcr-Chandrasekhar (dashed). The hatched regions 
are ruled out by requiring that less than 100% of WDs have com- 
panions with P<12 days (upper right corner) and that WDs do 
not merge faster than they are formed (bottom part of the plane) 
- see text for details. 

in galaxies of various Hubble types. For Sbc spirals with 
the stellar mass of the Milky Way (6.4 ± 0.6 x 10^° Mq, 
iMcMiUadlMll ) the specific SN la rate is 1.1 x 10"" SNe 
yr~^ Mq^- This number is remarkably close to our mea- 
sured specific WD merger rate, but an order of magni- 
tude higher than the super-Chandrasekhar merger rate. 

The implication is that there are not enough super- 
Chandrasekhar DD systems in our local region of the 
Milky Way to reproduce the measured SN la rate through 
the classic DD channel. Several authors have already 
pointed out the apparent dearth of super-Chandrasekhar 
SN la progenito rs usi ng qualitative argu ments (e.g., 
Ilsern et al.lll997t lMaog[2008l: iRuiter et al.l [20091) . Our 
analysis shows this from a quantitative measurement of 
pre-merger WD binaries. However, we find a remarkable 
agreement between the total WD merger rate and the SN 
la rate. For our assumed primary and secondary mass 
distributions, ~90% of mergers are sub-Chandrasekhar, 
but the total masses are often relatively high - 10- 
30% of the mergers are > 1.2 M©, and 25-50% are 
> 1.1 M0. Interestingly, recent theoretical work has ex- 
plored the possibility of su b-Chan drasekh ar CO /CO WD 
mer gers leading to SN la (van Kerkwiik et al.ll20l"ol see 



also iGuillochon et all l20ld iPakmor et al.ll20'Tl . Apart 
from the consistency between SN la rates and total 
WD merger rates, sub-Chandrasekhar explosions may 
have the advant age of producin g the correct chemi- 
cal stratification (jSim et al.ll2010f ). without r esorting to 
the ad hoc delayed detonation mechanism (jKhokhlovl 
|1991| ) needed by super-Chandrasekhar models. The mea- 
sured WD merger rate is also important for e stimatimg 
'foreg rounds' for gravitational wave detectors (jNelemanj 
[200l . 

Unavoidably, our conclusions depend to some degree 
on the assumptions made in the analysis. First, all WD 
samples are flux limited, which could lead to selection 
effects and unrepresentative values of /bin, et, or the WD 
mass distribution. For example, a DD population with 
a high super-Chandrasekhar merger rate may exist, but 
remain unobserved. Even so, the DD population we ob- 



serve, with its high merger rate, does exist, regardless of 
flux limits. Second, we have modeled the initial WD sep- 
aration distribution with a power law that is independent 
of mass. The true post-common-envelope distribution 
could have some preferred scale correlated with the com- 
ponent masses, which would affect the parameters and 
merger rates deduced from the ARVmax distribution. We 
have chos en primary masses b ased on the distribution ob- 
served bv lKeoler etall (12001 for SDSS WDs, which are 
mostly single, and we have drawn secondary masses from 
a flat mass-ratio distribution. These distributions might 
be different in DD systems. Different mass distributions 
would only affect weakly the binary population param- 
eters and the total merger rate, but they would have a 
strong effect on the super-Chandrasekhar fraction. We 
might have underestimated ARVmax for some double- 
lined WD binaries that are unresolved at the SDSS spec- 
tral resolution, and this would lead to an underestimated 
merger rate. Conversely, some of the high ARVmax sys- 
tems could conceivably be due to faint M star, rather 
than WD, companions. However, an M star of low 
enough mass to go undetected in the SDSS photome- 
try, < 0.15 Mq, and with a period of > 1 hour, would 
only induce a small ARVmax h^ the primary WD, e.g., 
220 or 270 km s-\ for WDs of 0.6 Mq and 0.4 Mq, re- 
spectively. Such systems therefore cannot dominate the 
high-ARVmax tail. At the spectral resolution of SDSS, 
the combined absorption plus emission line profiles of 
some WD+M binaries could mimic a large ARVmax, but 
this should not occur identically for several Balmer lines, 
so these systems would be removed by our vetting proce- 
dure. Finally, we have assumed that the velocity errors 
are normally distributed, with standard deviations quan- 
tified by our empirical tests, but there might be some low, 
non-Gaussian, tails to the error distribution that con- 
taminate the ARVmax distribution with false positives. 
In upcoming work, we will address these caveats by in- 
vestigating the effects of different input assumptions on 
the conclusions. We will also improve the purity of our 
ARVmax distribution with spectroscopic follow-up of a 
large number of binary candidates, and use this improved 
distribution to refine our constraints on the Galactic DD 
population and its merger rate. 
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TABLE 1 

Local WD Merger Rates and 95% Confidence Limits 







Total rate 


Super-Ch rate 


a 


/bin 


(10"^-' mergers yr^i Mq-^) 


(10^^^ mergersyr"^ Mq-') 


entire range 


0.014 to 0.32 


1.4 (0.16, 7.2) 


0.1 (0.016, 0.4) 


1.0 


0.11 to 0.24 


0.3 (0.065, 0.5) 


0.03 (0.017, 0.045) 


0.0 


0.046 to 0.22 


1.0 (0.46, 2.2) 


0.08 (0.03, 0.16) 


-1.0 


0.021 to 0.11 


3.0 (1.0, 6.0) 


0.16 (0.05, 0.3) 
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